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ABSTRACT 

The special role of the isoscalar mesons (a and us) in the NJL soliton is discussed. 
Stable soliton solutions are obtained when the most general ansatz compatible with van- 
ishing grand spin is assumed. These solutions are compared to soliton solutions of a 
purely pseudoscalar Skyrme type model which is related to the NJL model by a gradient 
expansion and the limit of infinitely heavy (axial-) vector mesons. 
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Introduction 

The theory of strong interactions, Quantum Chromo Dynamics (QCD), reduces for 
a large number of colors iV c to an effective theory of weakly interacting meson fields]]]]. 
Furthermore baryons emerge as soliton solutions in this effective meson theory 0. On the 
other hand, simplified models with quark degrees of freedom can also possess solitonic 
solutions, as for example the Nambu-Jona-Lasinio (NJL) model||. Selfconsistent soliton 
solutions have been foundH-0 where the calculations had been restricted to the chiral 
field. It has then been observed that a straightforward extension of this model to allow 
for space dependent scalar fields leads to a collapse of the soliton^ ||. Simultaneously 
the NJL model containing the chiral as well as (axial-) vector meson fields was shown to 
possess stable soliton solutions [pj]-]|l3|. In case the scalar degrees of freedom are incorpo- 



rated in the latter type of model, the repulsive character of the uo meson is supposed to 
prevent the soliton from collapsing. 

The aim of the present investigations is twofold. First, we will discuss the special role 
of the isoscalar meson fields related to the above mentioned collapse. With all meson 
fields included in the soliton calculations the valence quarks are strongly bound and they 
indeed join the Dirac sea. Therefore the topological current, which arises in leading order 
in the gradient expansion of the vacuum part of the baryon current, becomes a suitable 
representation of the full baryon current. This fact supports the Skyrmion picture of 
baryons. For these investigations we are going to employ a static energy functional 
which has recently been motivated by studying the relevant analytic properties under 
Wick rotation||14|. For the inclusion of all mesons similar calculations have previously 



been performed using a somewhat different definition of the energy functional [|Hi|j. As the 
Minkowski energy functional is ambiguous if the u meson is treated non-perturbativelyfln]] 
it is, of course, interesting to study the question of stability assuming an alternative 
definition for the energy functional. Combining the results of the two approaches should 
provide a general picture of the situation. 

On the other hand it is known that a Skyrme type model can be derived from the 



bosonized NJL model[|16|. This suggests to compare the self-consistent soliton solution 
of the NJL model to the soliton solution of the corresponding Skyrme model. Such a 
comparison constitutes the second issue of this letter. 

Description of the model 

As starting point we assume the bosonized version of the two-flavor NJL model IT 



A = A F + A m , (1) 

A F = Trlog(z J p) = Trlog(^+f+ 75 4-(P i ,S + P L St)), (2) 

Am = |rf 4 xf-^tr(StS-m (E + St)+m2) + i Ltr(^^ + A M ^)V (3) 

Here V^ = Y,l=o V^ ra /2 and A^ = Y,l=o A a ^T a /2 denote the vector and axial vector 
meson fields. The matrices r a /2 denote the generators of the flavor group (r° = 1). 
Pr,l = (1 ±75)/2 are the projectors on right- and left-handed quark fields, respectively. 
The complex field S describes the scalar and pseudoscalar meson fields S!y = X^ =0 S a T^j1 

E = S + iP = $ U (4) 

where we have made use of a polar decomposition. This defines the chiral radius $ as 
well as the chiral field U = exp(iO) with O being the chiral angle. The current quark 



mass matrix m — diag(mQ,mQ) only appears in the mesonic part of the action, A m . 
The Schwinger-Dyson (or gap) equation relates fn® to the constituent quark mass matrix 
rh = diag(m", m d ). For the ongoing discussion we will adopt the isospin limit m^ = 
itlq =: irtQ which also implies m u = m d = : m. As the NJL model is not renormalizable 
it needs regularization. This introduces one more parameter, the cut-off A. Hence the 
model contains four parameters: mo, <7i,<72 an d A. As ingredients from the meson sector 
we use the pion decay constant f n = 93MeV and the masses of the pion and the p-meson, 
m n = 135MeV and m p = 770MeV, respectively. Then the Schwinger-Dyson equation 
allows one to choose the constituent quark mass m as the only free parameter and express 



rn 0i 9i, 92 and A in terms of it. We employ Schwinger's proper time description [T^] which 
has the desired feature of preserving gauge symmetry. Then the fermion determinant, 
Af, m ay be decomposed into real (Ar) and imaginary (Ai) parts. Only Ar is divergent. 
The proper time procedure is applicable to Ar since PePe is positive definite 



Ar = ^Trlog (pip E ) - -1 £ ^exp {-sp^) . (5) 



1„ . /^ + ^ \ 1 f°° ds 

'l/A 2 

Although the imaginary part is UV finite, the proper time regularization may be applied 



as well|T2|, |13| . Note, however, that regularizing Aj or not yields qualitatively different 
models and different results may occur. 

The definition of the Euclidean Dirac one-particle Hamiltonian h via 

ippE = -d T -h (6) 

is useful in the context of studying soliton solutions since one has [d T ,h] = for static 
fields. In the limit of large Euclidean times, T — ► oo, the temporal part of the trace 
may be carried out straightforwardly by calculating Gaussian integrals. To evaluate the 
trace over the remaining degrees of freedom we require the eigenvalues of h. As h is 
non-Hermitean, we have to distinguish between left and right eigenstates 

h\m v ) = e v \m u ) (y v \h = e„(* v \ i.e. tf\* v ) = e* v \V v ) (7) 

with the normalization condition (^l^) = 5^. The functional trace may then be 
expressed as sums over the one-particle energies e M . The resulting expression motivates 
the following definition for the Minkowski energy functional [14]]: 

E[<p] = E^ + E^ + E* c + Ei ac + E m - E^ vac ] (8) 

where for simplicity we have generically labeled the meson fields by ip. The contribution 
of the trivial vacuum has been subtracted whereby the configuration <£> vac corresponds to 
E = m while the (axial-) vector meson fields are set to zero. E[tp] receives contributions 
from the explicit occupation of the valence quark level [T9|, EH] 

<i = N c J2vA^\, K,\ = N c J2VusMe*)el, (9) 

V V 

the polarized Dirac sea 

EL = ^El£|r( - \ (^/A) 2 ), EL = Z^L^elsgatf) (10) 

and the mesonic part, E m , which stems from A m - In eqn. (|) r] u denote the occupation 
numbers of the valence quark orbits. They are adjusted to provide unit baryon number: 



1 = J2n [Vfj. ~ 2 s S n ( e ^))- I n ec P- (03) we have displayed the vacuum part of the energy 
for the case that the imaginary part is not regularized; employing the regularization for 
the imaginary part yields |12|, |13| 



E L = ^E4sgn(ef)r(l (e*/A) 2 ). (11) 

It is important to note that the different treatments of real and imaginary parts under 
regularization destroys the analytical structure of the action in the time components of 
the (axial-) vector meson fieldsQ. As already indicated this prohibits a unique extraction 
of a Minkowski energy functional. In ref. |14]], however, it has been demonstrated that the 



definition (|8|) is motivated from the regularized energy functional in Euclidean space. This 
motivation is based on the fact that the energy eigenvalue, e^, may well be approximated 
by a sum consisting of two functionals of the meson fields. One of these two functionals is 
almost independent of u while the other has only a linear dependence on u, see tables 2.1 
and 2.2 of ref. fnfl . Furthermore, the energy functional (H) possesses the correct behavior 



under global flavor singlet transformations and it also yields the current field identities. 
A unique Minkowski energy functional can only be obtained when expanding the fermion 
determinant in terms of the time components of the (axial-) vector meson fields while 
treating all other fields as non-perturbative background fields (cf. ref. flT3|). 

Let us next construct the most general static Euclidean Dirac Hamiltonian in the 
grand spin zero sector. The grand spin operator is defined as the sum G = I + cr/2 + r/2 
with I being the orbital angular momentum, cr/2 the spin and r/2 the isospin operators. 
For the chiral field the well-known hedgehog ansatz 

U(r) = exp(ir • f6(r)) (12) 

satisfies the condition of vanishing grand spin while for the scalar field this can only be 
accommodated by a radial function 

$(r) =m(f)(r). (13) 

For the (axial-) vector meson fields we impose the grand spin symmetric ansatze: 

V° = w(r)V, V A a =0, V? = e ak *r k G(r), 

A° = 0, Al =0, A1 = f i r a F(r) +5 ia H{r), { ' 

where the indices a, i and k take the values 1, 2 and 3. Then the Euclidean Dirac 
Hamiltonian reads 

h = ex. • p + iu(r) + m0(r)/3(cos0(r) + 275T • fsinB(r)) 

+-(a x r)-rG(r) + ~(a-r)(T-r)F(r) + -{a ■ r)H{r). (15) 

With the ansatze ( |T2"D - ( |T^ ) the mesonic part of the energy is given by 

'm?fl 

TV J TV 



E 



4tt /rfrr 2 (^^[m(0 2 (r) - 1) + 2m o 0(r)(l - cos0(r))] 

+ (!^) 2 [ G 2 (r) + l F 2 (r) + F ( r )#(r) + lH\r) - 2cu 2 (r)}). (16) 
9v 2 2 ' 



a This is manifest in eqns. (|9|jlC|) by the explicit appearance of real and imaginary parts of the energy 
eigenvalues e^. 



-| Icy 

Here gy is the universal vector coupling constant, gy = (-^T(Q, tjt)) which is related 
to the coupling constant g-i via the p-meson mass g v = Ag 2 m 2 fllf . 

The equations of motion for the meson fields are derived by extremizing the static 
energy functional (|8|): 6E[(p]/S(p = 0. In addition to the equations of motion listed in 



ref. 13 we obtain for the scalar field 



<p(r) = — cos©(r) TT7rt r / -; — (cos0(r) + 275T • r sin 6 (r))p(r, r). (17) 

m raiii J An v ' 

with the scalar density matrix p(a?, y) (which is bilinear in the eigenfunctions of h) being 
defined in ref. Jl3| . The formal structure of the other equations of motion is not effected 
by the presence of the scalar field 0. 

In order to diagonalize h and to solve the equations of motion 5E[ip]/5ip = we 
discretize the Hamiltonian (|i5|) in a suitable basis using a spherical cavity of finite radius 
D. Typical values are D = 4 . . . 6fm. For details about the numerical method see ref. [pl[ 



The numerical solution to 5E/5ip = on the whole range < r < D is plagued by 
finite size effects [[L3|]. In order to avoid these effects we substitute the exact solution to 



5E/5(f = by the phenomenologically motivated large distance behavior of the meson 
fields for r > r m (D/A < r m < D/2). For the chiral angle this can be extracted from 
the solution of the free (P-wave) Klein-Gordon equation: 

e( r ) r ^ 0(r) ~ L + -. (18) 

r V r J 

For the other fields we make use of their relation to the chiral angle in a local approxi- 
mation yielding the asymptotic behavioif] 

G(r) ~ 2 (r), 

w(r) ~ 0'(r)0 2 (r), 

<f>{r)-l ~ 2 (r) + ae- 2mr , 

F(r),H(r) ~ 0'(r). (19) 



The constants of proportionality, which are omitted in eqns. (|T8||T9| ), are fixed by making 
contact with the exact solution to 5E/5(p = at r = r m . For the scalar field we fix the 
additional constant a by requiring the derivative of <p(r) to be continuous. 

The role of the isoscalar meson fields 

In some recent studies of the NJL model J7|, || it has been shown that abandoning the 
chiral circle condition S 2 + P 2 = m 2 leads to the collapse of the soliton configuration if no 
additional mechanism to preserve the stability is incorporated (e.g. including a forth-order 



term in the scalar meson field |22|, |23[ or constraining the regularized baryon number 
It has already been demonstrated that in an alternative definition of the (ambiguous) 
non-perturbative energy functional the incorporation of the u meson provides a further 
stabilisation mechanism [HJ]. This is intuitively clear because the uj meson is of repulsive 
character. 

As already remarked in the preceding section regularizing the imaginary part of the 
action or not leads to quite different models. Taking a finite cut-off for the imaginary part 



b In order to accommodate the vacuum values of the fields at r = D wc furthermore multiply appro- 
priate factors to the RHS of eqns. (118) and (|l9|): ip — * <^tanh(a(l — j^)) with awlO- 15. 



Table 1: The soliton energy E as well as its Dirac sea and mesonic contributions E vac 
and E m for different values of the constituent quark mass m. Also shown is the energy 
of the 'dived' level e va i and the axial charge g A - 



m (MeV) 


350 


400 


500 


E (MeV) 


1125 


1091 


1022 


Etc (MeV) 


2271 


1337 


883 


EL (MeV) 


177 


206 


223 


E m (MeV) 


-1323 


-451 


-83 


e val/ m 


-0.28 


-0.51 


-0.72 


4al/ m 


0.15 


0.15 


0.12 


9a 


0.41 


0.36 


0.32 



TT|) also yields a regularized baryon number. As a matter of fact this regularized baryon 
number density vanishes for the field configuration associated with the above mentioned 
collapse [Q] indicating that the collapse is due to the transition from the baryon number 
Nb = 1 sector to the Nb = sector. Since the baryon number density represents the 
source of the u field nothing prevents the u field from being zero when the imaginary 
part of the action is regularized and the scalar field is allowed to be space dependent. 
Consequently, the collapse also appears if all vector mesons are included and Ai is reg- 
ularized. The ongoing explorations will therefore be constrained to the case when Ai is 
not regularized. Then a non-vanishing source for the u field is present and the stability of 
the soliton depends on the strength of the u field. It has been noted that the numerical 
calculations will pretend to a "pseudo-stable solution" if too small a basis for diagonaliz- 
ing h in momentum space is adopted ||. By varying the size of this basis we have ensured 
that our solutions are not subject to this "pseudo-stability". At low constituent quark 
masses, m, the trivial minimum with N c free and unbound valence quark orbits occupied 
is energetically favored compared to the soliton configuration. Actually a lower bound 
exists for m below which the equation of motion are solved by this trivial configuration 
only. Numerically we obtain for this bound approximately 330MeV which it is somewhat 
larger (~ 410MeV) when neither p nor ai fields are present. It should be stressed that 
this instability is completely different in character from the above mentioned collapse. In 
case of the collapse the valence quark is strongly bound and its energy eigenvalue tends 
to —m. 

For stable solitons the valence quark is strongly bound (e^j ^ — m/4) and the chiral 
radius <fi lies near the chiral circle <f>(r) = 1 (c/. fig. 1). Then the scalar-isoscalar meson 
field has only minor influence on the energy and on the other meson fields. In table 1 the 
soliton energy E, its Dirac sea and mesonic contributions, _E vac and E m , as well as the 
energy of the valence quark level (e val ) are shown for different values of the constituent 
quark mass m. Also displayed is the axial charge g A obtained directly from the axial 

J drr 2 [H(r) + F(r)/3]. The relative contributions from 



vector field 113 



g A = -(2ir/g 2 ) 

the Dirac sea and the mesonic part of the scalar-isoscalar meson field depends on the choice 
of the point r m where the tail for the meson field is fitted. Nevertheless the total energy 
and the valence quark energy of the soliton are stable in a wide range D/A < r m < D/2. 



To be definite, we choose r m such that d r 



which yields r m « 0.35 ■ • ■ 0.375D. 



Comparison with the Skyrme model 

Since the energy of the valence quark orbit is negative the baryon number can well 



be approximated by the topological current. For that reason the NJL soliton strongly 
supports Witten's conjecture that baryons may be described as solitons within purely 
mesonic models. We will now study the explicit connection to the Skyrme model. In a 
first step we assume a gradient type expansion for the pseudoscalar field. Secondly, we 
adopt the static limit for all other mesons, i. e. the inverse propagators are approximated 
by the corresponding mass termsB This allows one to integrate out these mesons yielding 
an extended Skyrme lagrangianFj 

£ = £2 + £4 + £-6 + £sb, 



£ 



£2 
£1 

£ e 

SB 



f 2 

|gtr [L M , L v ] [£/*, LI + S (tr {L^)f 

2 ^ 
flint 



24e 2 



tr 



(U + U ] - 2) . 



(20) 



Here L^ = U^d^U and B^ = (l/2A-K 2 )t^ up xtiL v L p L x denote the left Maurer-Cartan form 
and the topological current, respectivly. In the static limit the parameters e, e^ and x 
are obtained to be 



\2a 



Cfi 



tanx 



N 

N c g v 
2m p 



g v a 



+ 



V3A 

Am 



with 



1 H -, 

6m 2 ' 



g-ya 



4m |2a-l| 



(21) 



For a given chiral field L the soliton energy or mass obtained from the Lagrangian ((20|) 
is given by 



M 



d A r 



fl (o*,^ , 2sin 2 (9(r))\ 1 sin 2 (9(r)) / a ,. , sin 2 (6(r))' 



9"(r) + 



0"(r) + 



2r 2 



cosx 



-_IK) + ^ff 



smx 



'6 cv/2 



57T 



9"r 



sin 4 (9(r)) 



+ m 2 / 2 (l-cos(e(r)) 



.(22) 



If we perform slow rotations A(t) e SU(2) on the static soliton [/(£) = A(t)LMJ(t) we 
can derive the following Hamiltonian in the space of collective angular degrees of freedom 

J 2 



H = M + 
with A being the moment of inertia given by 



2A 



(23) 



- / tfVsin 2 (0(r)) 



f 

■J n 



2 , COS * Ufl 



6 /2 M 



sin 2 (6(r))' 



2sin% 
3e 2 



9 ,2 (r) + 



2sin 2 (0(r))\ 



+ 



4vr 4 



6 e /2 (r 



sin 2 (9(r)) 



(24) 



c In the NJL model the mass of the scalar meson is given by 2m. 
d For notation see ref. [E5| 



Table 2: Comparison of the mass M, eq. (p2|), the masses of the nucleoli (N) and A- 
resonance (^) for the self-consistent NJL model and the soliton of the extended Skyrme 
model (pCj). Shown are the results for two different values of the constituent quark mass 



m. 



The physical p meson mass is adopted 



, m p 



m 



ph 



m (MeV) 


400 


500 


soliton 


eq. © 


NJL 


eq. O NJL 


e 


4.65 


5.41 


e 6 (10~ 3 /MeV) 


12.7 


13.7 


X 


0.221 


0.191 


M (MeV) 


1760 


3519 


1714 


4124 


M N (MeV) 


1812 


3643 


1770 


4261 


M A (MeV) 


2023 


4139 


1995 


4809 


M A - M N (MeV) 


210 


496 


225 


548 



The mass for the N and A might be extracted from eq. (|23|): 

3 



M N 
M A 



M + 
8A 



(25) 



For comparison we have calculated the chiral angle 0(r) (cf. fig. 2): on one hand the self- 
consistent NJL soliton and on the other the extended Skyrmion obtained from the Euler- 
Lagrange equation associated with the Skyrme Lagrangian (p0|). We have then calculated 
the N and A masses (E5f) for the NJL soliton chiral field as well as for the Skyrmion. As 



can be seen from fig. 2 using the physical p meson mass m p = m ph = 770MeV leads to 
quite different profiles for the chiral angle 0(r). Accordingly, the soliton mass M and 
therefore the N and A masses are quite different in the static approximation(c/. tab. 2). 
Note that the soliton mass M is significantly larger than the self-consistent NJL soliton 
energy (see table 1). This indicates that the physical p meson mass is too low for the 
validity of the static limit. Indeed if we increase the p meson mass successively the two 
chiral angles approach piece by piece (cf. fig. 2). For m p = 8m ph the soliton mass M 
obtained by substituting the self-consistent chiral angle of the NJL soliton into the static 
approximation and the Skyrmion are in very good agreement as can be seen from tab. 3. 
Furthermore, this soliton mass is very close to the corresponding self-consistent NJL 
soliton energy which is E = 1162MeV (E = 1026MeV) for m p = 8m p p h and m = 400MeV 
(m = 500MeV), respectively. 



Conclusions 

To summarize, we have shown that for an appropriate definition of the NJL model 
energy functional the inclusion of all meson fields in a way compatible with the hedgehog 
ansatz leads to stable soliton solutions in a wide range of constituent quark masses. We 
would like to put emphasis on the fact that the chiral radius deviates only mildly from 
its vacuum expectation value when all meson fields are included in the evaluation of the 
self-consistent soliton. This justifies the frequently adopted approximation <p(r) = 1. 

Including all meson fields with vanishing grand spin the valence quarks are strongly 
bound and the baryon number is completely carried by the polarized Dirac sea. This 
indicates that baryons can be described as purely mesonic topological solitons like the 
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Table 3: Same as table 2 for m p = 8m^ h 



m (MeV) 


400 


500 


soliton 


eq. © 


NJL 


eq. (|0|) NJL 


e 


5.17 


6.46 


e 6 (10" 3 /MeV) 


2.53 


3.15 


X 


0.274 


0.274 


M (MeV) 


1160 


1162 


1008 


1013 


M N (MeV) 


1326 


1328 


1235 


1272 


M A (MeV) 


1989 


1991 


2140 


2307 


M A - M N (MeV) 


663 


663 


906 


1035 



skyrmion. This result has motivated the comparison between the solitons of the NJL 
model and a Skyrme type model. The latter has been related to the NJL model in the 
limit of large (axial-) vector meson masses and a gradient expansion for the pseudoscalar 
field. It has turned out that these two models possess different soliton solutions for 
the physical values of the vector meson mass (m ph ). However, as m p is increased these 
two solitons become similar in shape and size. Reasonable agreement is achieved for 
m p > 8m ph . This indicates that the kinetic term for the (axial-) vector mesons, here 
represented by the fermion determinant, carries important information about the meson 
fields. We have thus collected some support for Skyrme type models from a microscopic 
model for the quark flavor dynamics. One might argue that purely mesonic models with 
explicit (axial-) vector degrees of freedom might be more feasible for the computation 
of nucleon properties. However, as one is interested in the explicit quark structure of 
baryons it is unremitting to consider a microscopic model like the one of Nambu and 
Jona-Lasinio. 
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Figure captions 



Figure 1 

The chiral radius field d>(r) for m = 400MeV. 



Figure 2 

The chiral angle 0(r) of the NJL model and the extended Skyrme model for m = 400MeV. 
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